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A new diagrammatic method, which is a reformulation of Berezinskii's technique, is constructed 

to study the density of electronic states p{e, (f>) of a one-channel weakly disordered ring, threaded by 

^^ \ an external magnetic flux. The exact result obtained for the density of states shows an oscillation 

of p(e, 0) with a period of the flux quantum 00 = ^ ■ As the sample length (or the impurity 

^ ' concentration) is reduced, a transition takes place from the weak localization regime {L ^ I) to the 

^J( , ballistic one (L < I). The analytical expression for the density of states shows the exact dependence 

^^ ■ of p(e, 0) on the ring's circumference and on disorder strength for both regimes. 

cn 

71.10.-w, 71.20.-b, 71.55.i, 73.23.-b 



The oscillation of physical properties of disordered metals has been studied intensively after the prediction of the 
Aharonov-Bohm effect in doubly connected dirty systems W with the period of half of a flux quantum and its 
^ ■ observation [g[ in a Mg cylinder. 

Today, a particular subject of intensive investigation is the persistent current, predicted in pMl for one-dimensional 
^N . disordered rings. Recent advances in microstructure technology facilitate the fabrication of mesoscopic rings and the 
observation of thermodynamic currents therein |^-0|. The observed oscillatory responses in these experiments, which 
are consistent with a persistent current, differ in the period of oscillation. 

A similar controversy exists also in theory. According to fundamental physical principles all physical parameters, 

in particular the persistent current, of a one-channel metal ring should be periodic in an applied magnetic flux 

S ■ with period of a flux quantum (t>o — ^ |^,^,H-|0[. However, the coherent backscattering mechanism with consequent 

i-rt ' interference effects in mesoscopic systems gives rise to conductance oscillations with the halved period 0o/2 p|,pT|-p^. 

r^ 1 It is pertinent to notice that the attempt to explain the (t>o/2 oscillation in a disordered ring by taking into account 

O ' the electron-electron interaction p5|-pj| is also based on the "cooperon" propagation in the system. 

, ^, , All these disputes in the theory seem to be connected with the absence of a consistent theory for a one-dimensional 

(Id) disordered ring in a magnetic field which goes beyond the diffusion approximation and can calculate not only 

P\| ! average values of the physical parameters but also mesoscopic fluctuations of these parameters. 

^ ' It is well known that the physical parameters of a mesoscopic system with dimension L satisfying the condition 

C^^ ' I < L -^ lin (where I is the mean free path and ^in is the length over which the phase coherence of an electron wave 

\0 ' 

O 

O 



is conserved) have random character, i.e. self-averaging is violated |20|. At T = all systems become mesoscopic. In 
this case high moments give a considerable contribution, which results in strong differences between average value and 
typical one of the observed parameter |Q, i.e. the average value loses its significance to characterize the experimental 
qn^ ' observation. For such a problem one has to calculate the whole distribution function and to get the typical value for 
^^ an observable parameter [pi| , p2[ . 

^^ To perform the procedure presented above there exist technical difficulties. As far as the Aharonov-Bohm problem 

2 , for a sufficiently narrow ring is Id, the diffusion approach does not give correct results because of strong interference 
^ • effects independent of the degree of randomness [p3 . The periodicity adds an additional technical difficulty. 
' In this paper we present a new diagrammatic technique by means of which all diagrams can be summed exactly 

y ■ for weak disorder, when the criterion kpl ^ 1 (where kp is the Fermi momentum) is satisfied. This method is a 
r^ ' generalization of the Berezinskii method pj| , which was previously developed for a strictly Id system. 
^ , The latter system with J-correlated Gaussian impurity potential was studied a long time ago by Halperin ||25| . In 

L| ' difference to our case, Halperin considered the limit of an infinite density of scatterers, where the loffe-Regel criterion 
. ,—1 I {kpl ~ 1) is reached. Halperins result describes the energy dependence of the density of states (DoS) of bound states 
j^ ' appearing in the impurity tail with negative energies. The same results for the DoS, together with new information 
%-i ^ on the localization length, dielectric constant, and conductivity, were later obtained in an extension of Berezinskiis 
. . . , theory to strong disorder, [ p6| , p7[ . 

Here, we consider a one-channel metal ring, threaded by a constant magnetic flux 4> through the opening. The 
electrons inside the ring with circumference L are elastically scattered through the impurity potential Vimp(a;)- The 
Hamiltonian of the system is written in the form 

TT ^^ /■ d 2Tr (j) ,2 ^, , . ,,, 

where x — ip-^;^ is the spatial variable on the ring, <j)o = -^ is the fundamental period of a flux quantum and m* is 



the effective mass of an electron. The impurity potential l^mp {x) is considered here to be Gaussian distributed with 
a spatial width small enough to justify the Born approximation. We apply here our new diagrammatic method to 
study the DoS at T = according to the expression 

p(e>;a;) = --Ini(G+(x,x;e)} (2) 

TT 

where G'^{x,x';e) is the retarded Green's function (GF) and the bracket means averaging over the impurity realiza- 
tions. 

Berezinskii's idea to construct a real space diagrammatic method in one dimension is based on the factorable 
form of the "bare" GF. However, for a Id problem with periodic boundary condition, the quantization of the energy 
spectrum creates difficulties in this respect. To avoid these difficulties, the boundary condition is not imposed at the 
beginning and we start with a free particle of energy efe = T^i^ ~ X"^)^ ^'^'^ wave function ^^(a;) oc exp(ifca:) with 
continuous k. To implement the periodic boundary conditions, the particle is allowed to make an arbitrary number of 
revolutions around the ring in both directions. By this means, the "bare" retarded Green's function in the coordinate 
representation can be expressed in factorable form, as it takes place in the Berezinskii technique p4| 

J 2tt e-ek±iri nv{e) \ (po L v[e) J 

where the dissipative parameter rj characterizes an energy level broadening due to inelastic scattering, v{e) and 
tip{e) are the velocity and the momentum of an electron with energy e in a strictly Id system, respectively, with 

v{e) = \/^, pie) — \/^w^- It is worthwhile to note that the true GF for a clean ring, Gq , can be obtained from 



Eq.(0) by making allowance for arbitrary revolutions. Gq{x, x'; e, (f>) can then be expressed in terms of Gq as 

oo 

G+(x,a;';e,0)= ^ G+(a;, x' + riL; e,</.) (4) 

n— — OO 

We can easily verify according to Eqs.(||)-(U) that (the impurity averaging in Eq.(||) of course loses meaning in this 
case) the DoS of a clean ring in the presence of an external magnetic field is 

oo , 

po{e, (p) ^ po + 2po V" cos{p{e)Ln) cos(27r-— n)e"^^^" (5) 

where the DoS of a strictly Id system is denoted by po = TrvU)h ' ^^ ^^^ limit 77 ^ 0, the DoS assumes the expected 
discrete form. 

Now we represent the unaveraged retarded GF G^(x, x'; e, (f)) of an electron moving in a field of randomly distributed 
impurities by a continuous line going from point x to x' . For the DoS problem, G~^{x,x;e,4') is adequate. The 
factorable structure of the "bare" GF Gj(a;,a;'; e, (^) between two subsequent scattering points Xi and Xi+i makes 
it possible to transform the coordinate dependence from the line to the impurity scattering points Xi, where the 
impurities are located. For our case of weak disorder, the diagrams have to be selected taking into account ppl ^ 1, 
where pp is the Fermi momentum and I is the mean free path. The correlator, connecting two scattering events 
and depicted in the diagrams as a wavy line, characterizes the essential vertices shown in Fig.|l|. The expressions 
corresponding to these internal vertices are a) — [-^ + jjr)? b) ~it: ''^^'^ c) —-j^, respectively, where l'^ and l~ are 
the mean free paths with respect to forward and backward scattering [p7| . The contribution of the neglected vertices 
vanishes, since they oscillate strongly with the position |24]. These expressions show that the internal vertices don't 
depend on the magnetic flux (j) and the magnetic field dependence can be transfered from the internal vertices to the 
external ones, which are given in FigJ2. 





a) b) c) 

FIG. 1. The three internal vertices giving an essential contribution to the impurity averaged GF in the weak disorder limit 

ppl > 1. 
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FIG. 2. The external outgoing (a,d) and incoming (b,c) vertices and the expressions corresponding to them. 

In difference to two-point correlator problems, the DoS problem can be described by rather simple diagrams, an 
example of which is presented in Fig.|3|a. It is convenient to cut the diagrams at point x and straighten the lines to 
arrive at the form shown in Fig.0b. 

Each diagram is characterized by the number of line pairs m returning to the cutting point x and by the total 
number of throughgoing lines n = n+ + n~, where n+ and n~ are the numbers of rightgoing and leftgoing lines, 
respectively. According to Berezinskii's method ^4| , the sum of all diagrams having m pairs of returning lines and n 
throughgoing lines at the cross section x is denoted by Qo^m, n; x — x'). This block does not contain the contributions 
of the external incoming and outgoing vertices, which are included in the final expression as additional multipliers. 
Due to the structure of the essential vertices, the number of loops on the left hand side is identical to the number of 
loops on the right hand side and we can use the common symbol m. 




X X x+L 

FIG. 3. (a) A diagram giving a contribution to the DoS. The radial unfolding of the drawing was done for the sake of clarity, 
(b) The same diagram as in (a) after cutting at the point x. It belongs to the class of diagrams with (m = l,n+ = l,n~ — 0). 



The expression for the average value of the retarded GF can be written as 
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(6) 

where the combinatorical factor in the angular brackets denotes the different possibilities of ordering the loops and 
lines: assuming that the electron starts from the left hand side, and pursuing the continuity of an electron line for the 
GF, the 77+ rightgoing lines can be distributed arbitrarily on the 777 + 1 positions before each of the loops on the left 
side and directly before the final external vertex. Also, the 77~ leftgoing lines can be distributed on the 777 positions 



before the loops on the right side. This gives the first term in the angular brackets of Eq.(g) from 
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(7) 



Similarly, the electron can also start from the right hand side. In this case, the left- and rightgoing lines reverse their 
role, and this makes the second term in the angular brackets of (m). However, if there are no loops and no lines, these 
two cases can not be distinguished, therefore one has to include the (5~term in the angular brackets as compensation. 
The exponential term in Eq. (H) comes from the external vertices by taking into consideration the circulations around 
the ring. 

For the final expression for the DoS, we insert Eq.(pf) into Eq.(0). After some transformation of variables we obtain 
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(8) 



The equation for the central block Qa{m, n; x) is constructed by infinitesimal shifting the point x and examining the 
change in Qo due to passing of the individual impurity lines through x [E4U2j] . This process is schematically presented 
in Fig.Q. The numbers of possible insertions of the vertices a) and b) of Fig]l]are (27n, + 77) and ^{2m + n){2m + n— 1), 
respectively. The vertex c), however, can be inserted in two different ways: i) without changing 777 and 77, this can be 
done in 777(777 + 77 — 1) different ways (First two blocks in Figfl); and ii) with changing 777 to 777 — 1 and 77 to 71 + 2 
(Third block in Fig.0). The latter way of insertion has 777^ possibilities. 
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FIG. 4. Scheme to construct the equation for the central block Qo{m,n; x) 
In total, the equation for Qq is 



dx 



Qo{m,n;x) = 

/ (2777 + 77)^ 77 777(777 + 77) 
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Qoim,n;x) 



—777^ Qo{m — 1, 77 + 2;2;) 



Qo satisfies the boundary condition 



)o(777,77;a; = 0) = djn,o 



which means the absence of scattering for a ring with an infinitesimal small circumference. 
To solve Eq.(0), we replace QQ{m,n;x) according to 



Qoirn, 77; x) = exp(— — — (2777 + n)^ — —m{m + n) — -;rf^n)QQ{m, n; x) 



21 



(9) 



(10) 
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By Laplace transforming Qo from the coordinate x to the new variable A and by using the boundary condition (fKl), 
the equation for Qo is reduced to the form 



1 'VJl Tl 

XQQ{m, 77; A) - (5„i^o = --r^rri^ exp(— )(5q(777 - 1, 77 + 2; A - — ) 



(12) 



Eq.(p^ can be solved by iteration in m. For m = 0, Qq{0, n; X) — j. Further iteration gives 

Qo{rn,n;X) - ^ ' _ "^ 11 x i ■/ ^ ^ (13) 

Inverse Laplace transform of Qf^{m, n; A) results in 

Qo{m, n; L) = exp(- — (27n + n^ - —m[m + n) - T^n) 

m\ m\{j + n-2)\ ■, ( ^ - (^^^ 



j=o 






where the exponential prefactor in ( pi] ) has been taken into account. From Eq.(p^ it can be verified that for L = 
the sum over j gives Qo{m,n, L = 0) = Sm,Q and that Qq decays exponentially with n for m = Q. Also, one obtains 
from (pi) in the special case of 7i = 



)o(m,n = 0;L) = (1 -m-m— )e TT i- + e TT ,- «- N (-l)J — [2j - l)eii-'- ■' ' 

I- ^ \] [m+j -1)\ 



(15) 



Eqs.(g|) and (|14|) constitute the exact result for the DoS of a one-channel weakly disordered ring in an external 
magnetic field. The result is valid for weak localization and ballistic regimes. 

For the weak localization regime, corresponding to the criterion L 3> max{/+,^^}, Eqs.(^ and ( |l4| ) are simplified 
to 

p{eA) = Po{l - p exp(-p - A _ _!L_i)} 

+2pocxp(- A _ _L _ _!L_L){cos(p(e)L)cos(27r^) +exp(-|^ - A) cos(2p(.)i) cos(4.|-)} 

which shows that the leading contribution to the DoS oscillation has a period of 0o and its amplitude decreases expo- 
nentially with impurity strength (or with increasing L) for weak disorder. Such a small contribution of the impurity 
scattering to the DoS is connected with the absence of "diffusion" and "cooperon" contributions to the averaged 
Green's functions. In the limit of an infinite sample, the correction due to weak disorder disappears completely. For 
the ballistic regime, when L < imn{l~^ , Z~}, the contribution to the DoS can be approximated in the form 

p{e, (j)) = po{e, 4>) - po— V" n^ cos(p(e)Ln) cos(27r— n)e~^^-^" 

n=0 ^^ 

L ^' 6 r, 

— Po-;—}^ ncos{p{e)Ln)cos{2Tr-—n)e^'<^ " (17) 

/ ^-^ (hn 

- 2po73 V V(fc + 1) cos(p(e)Ln) cos(27r— (2/c - n))e"^^" 

where pQ{t,cj)) is the DoS of a clean ring given by Eq.(||), iV+ « [\/ Tt] ^^^ ^- ~ [tt] • From Eq. ([l7|) it can be 
seen that disorder gives a contribution proportional to j^ to /Oo(e, 0). the DoS of a clean system in a magnetic field, 
Po{e, (f) — p() + 2po X]^i cos(p(e)in) cos(27r-^ri). The dependence of the DOS on the ring circumference L, computed 
from Eqs.(|^) and (|l4|), is given in Fig.g|, where the sharp discrete peaks in the ballistic regime are due to energy level 
quantization. 
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FIG. 5. Dependence of the DOS on the ring length L for smaU rings and zero magnetic field, obtained from Eqs.dsl) and 
([L4[). The fermi momentum is given by pl~ = 100. The sharp discrete levels [Eq.(|l7[)] for small ring length or large scattering 
length cross over to a continuous DOS for large rings [Eg. jig)]. 

The upper limit N± of the sums in Eq.(p^) may be a small value, e.g. N± « 1, deduced from -^ = 1.3 according to 
the experiment in Ref. B. Therefore, the oscillation with a full flux quantum 0o will be pronounced in the ballistic 
regime. 

In the absence of backward scattering (?^ — oo) in the system, Eqs.(Bh and (14) give a rather simple expression for 
the DoS, which can be presented in the following form: 



P(e» == Po + 
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exp[-ip(e)L - i27r(?!)/0o + ^7 + ih)^ - 1 exp[-ip(e)L + i2-K4)/(l)t) + 17 + 4^L] - 1 



c.c. 



(e)^ 



(18) 



Eq.(n8h can be physically interpreted as follows: each act of forward scattering gives rise to coherent shifting of all 
energy levels. The value of this shifting is random with Gaussian distributions; the typical value of this shifting is 

proportional to -^y -r where t~^ is the relaxation time due to forward scattering. Therefore, only backward scattering 
seems to be responsible for level repulsion in disordered Id systems pq |. Averaging over this random shifting results 
in Gaussian broadening of the energy levels, which for the weak localization regime is much smaller than Dingle 
broadening and comparable with it in the ballistic regime. In the latter case, the transport seems to be connected 
with resonant tunneling. 

To illustrate the dependence of p on Z+ and l^ , we decomposed the DOS [Eq.(||)] into the field independent part 
with the restriction An = |n+ — n^ | =0, and the harmonics with An — 1,2.... For /+ = /~ and for /+ = 00 we show 
these contributions in Fig.pl as functions of the ring length. For l~ — 00, the field independent part is constant and 
the higher harmonics are simple damped oscillations. 
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FIG. 6. The field independent (An = 0) part of the DOS for 1+ = r (l.h.s) and for /+ = cx) (r.h.s); and 
insets show the An = 1 contributions. Higher oscillations (An > 1) are similar, but with increased damping. 




50. The 



It is necessary to notice that the oscillative behavior of a persistent current will differ from that obtained for the 
DoS. In contrary to the dynamical approach to the study of conductance, as a result of which the latter is connected 
with a current-current correlator, the average thermodynamic current {I{(f>)) is defined by the average value of the 
thermodynamic potential F according to (/('/')) — — c i. where the bracket denotes an averaging over the impurity 
realizations. 

At zero temperature the last expression turns to (/(</>)} — ~c-4^ with E = J^ deep{e, cj)) being the total energy 
of the particles. Expressing DoS p(e, (p) and the flux dependent Fermi energy /x((/)) as p = (p) + Sp and p — (p) + dp 

where (Sp) — (Sp) = and using in addition the particle number conservation N = const = L dep{£,(l)) to 
determine Sp, (E) can be written in the following form: 

(E) = r de 6(p(e, 0)) - ^ r de {5p{e, <P)5p{p^,<f^)) 

Jo \P) Jo QqN 
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dti I de2{Sp{ei,(j))Sp{e2,(l))) 



As it is seen from this expression, contributions to the persistent current are given not only by the average value of the 
DoS but also by the correlator {6pSp). By expressing the DoS as a difference of the retarded (G"'") and advanced (G~) 
Green's functions, the latter correlator is shown to be dominated by {G~^G~). Characterizing the retarded (advanced) 
Green's function by the number of line pairs m (jn) and by the total number of throughgoing lines n — n^ + ?i^ 
(n = n+ +n~), the oscillating factors in the kernel for the correlator {G^G~) will have the form (compare with 
Eq.(p|)) cos(p(e)L(n — n)) cos(27r^(n+ — n^ + 71+ — rT))- The main contribution to the correlator comes from the 
terms with n = n, which oscillate with the halved period, p9| . 

The explanation presented above for the ^ ~ periodici ty of { I{4>)) is connected with the transition from the grand 
canonical ensemble averaging to the canonical one [p yil| , p^ , [l7| , [l8| , p0| . Another consideration also seems to be pertinent. 
For a mesoscopic system the total energy _E at T = is a random parameter. For a fixed chemical potential in the 
system the number of levels fluctuates from sample to sample. An external magnetic field will periodically change 
these fluctuations. The correct approach to the problem is to find the distribution function of E, which seems to be 
the same as the one for the DoS, and to calculate the typical value of E by means of averaging. Such averaging will 
include not only (p) but also higher moment correlators. 

The method presented in this paper makes it possible to calculate all these moments of the DoS [M, as it was 
done in |2^ for a strictly Id system. Further, the extention to correlators of the local DOS with different energies 
and different positions allows to study level repulsion. 

In conclusion, we have presented a new diagrammatic method which gives an exact result for a one-channel weakly 
disordered ring threaded by a magnetic flux <f>. As an application of the method the DoS of an Aharonov-Bohm ring 
in the absence of electron-electron interaction is calculated. The result obtained gives an exact dependence of the 
DoS on the parameters 0, e, L and I . 

We are thankful for useful discussions with V. N. Prigodin in the early stage of this work. This work was supported 
by the SFB410. 
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